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Pulsars are rotating, magnetized neutron stars that are born in supernova explosions following the 
collapse of the cores of massive stars. If some of the explosion ejecta fail to escape, it may fall back onto 
the neutron star^ or it may possess sufficient angular momentum to form a disk.^ Such 'fallback' is both 
a general prediction of current supernova models'' and, if the material pushes the neutron star over its 
stability limit, a possible mode of black hole formation.^ Fallback disks could dramatically affect the 
early evolution of pulsars, ^'^ yet there are few observational constraints on whether significant fallback 
occurs or even the existence of such disks. Here we report the discovery of mid-infrared emission 
from a cool disk around an isolated young X-ray pulsar. The disk does not power the pulsar's X-ray 
emission but is passively illuminated by these X-rays. The estimated mass of the disk is of the order 
of 10 Earth masses, and its lifetime ( >, 10^ yr) significantly exceeds the spin-down age of the pulsar, 
supporting a supernova fallback origin. The disk resembles protoplanetary disks seen around ordinary 
young stars, ^ suggesting the possibility of planet formation around young neutron stars. 



The so-called anomalous X-ray pulsars^ (AXPs) are 
a group of young ( ^ 10^ yr) neutron stars with spin 
periods falhng in a narrow range (5-12 s), no evidence 
for binary companions, and whose X-ray luminosities 
(~ 10'^^ erg s~^) greatly exceed their rates of rotational 
kinetic energy loss (~ lO'^^ erg s~^). AXPs are gener- 
ally believed to be 'magnetars',^ which are isolated neu- 
tron stars with exceptionally strong ( ^ 10^"* G) surface 
magnetic field strengths and whose magnetic energy ulti- 
mately powers their X-ray emission. An alternative ex- 
planation for AXPs attributes their X-ray emission to ac- 
cretion from a residual debris disk,^"^^ but this model has 
had difficulties explaining observations.^^ The brightest 
known AXP is the 8.7 s pulsar 4U 0142-f61, at a distance 
of 3.9 kpc (ref. 13). Besides its X-ray emission, the pul- 
sar also has known optical^^ and near-infrared^^ (near-IR) 
counterparts. 

As part of a systematic search for fallback disks 
around young neutron stars, we observed the field around 
4U 0142-1-61 in the 4.5 /im and 8.0 fim bands with the 
Infrared Array Camera (IRAC) on the Spitzer Space Tele- 
scope to look for the infrared excess predicted by models 
for an X-ray heated fallback disk.^^ We found a candidate 
mid-IR counterpart at the pulsar's position in both bands 
(Fig. 1) that has very unusual IR colours (Fig. 2). Based 
on the position coincidence and colours, we conclude that 
we have identified the mid-IR counterpart of 4U 0142+61. 

We can reconstruct the observed low-energy spectral 
energy distribution of 4U 0142+61 by combining our 
Spitzer data with existing IR and optical data^^'-'^^ (Fig- 
3). We may then infer the intrinsic spectrum by correct- 
ing for interstellar reddening. Although this reddening 



correction has little effect on the mid-IR data, it signifi- 
cantly affects the optical data. For any reasonable choice 
of reddening (2.6 < Ay < 5.1; ref. 12), the optical and IR 
data clearly arise from two different spectral components. 
The optical (VRI) emission is consistent with a power-law 
spectrum, ^^ is pulsed at the spin period,^'' and is presum- 
ably of magnetospheric origin; we will not consider this 
component further but will instead confine ourselves to 
the distinct IR component. Hereafter, we adopt the most 
likely reddening value of Ay — 3.5 (ref. 13). 

We begin by noting that a non-thermal (self-absorbed 
synchrotron) origin for the IR component can be ruled 
out, because associating the turnover near 4.5 fim with 
a synchrotron self-absorption break requires an emission 
region with an implausibly weak (~ 10^ G) magnetic field 
strength within a small (~ 100 km) radius of an active, 
highly magnetized pulsar. Similarly, the shape of the 
IR spectrum suggests a thermal origin, but the simplest 
case of a single-temperature (~ 920 K) blackbody gives 
a mediocre fit to the data. Moreover, the only plausible 
blackbody source would be a large planet or extremely 
low-mass stellar companion, but the implied blackbody 
emission radius (> 5Rq, where Rq is the solar radius) is 
much too large. 

A multi-temperature (700-1200 K) thermal model fits 
the data better and thus naturally suggests an extended 
disk or shell origin for the IR emission. We reject a shell 
geometry as it would produce ^ 1000 times the X-ray and 
optical absorption observed toward the AXP. As for a disk, 
it has already been shown that an accretion disk model 
(where the pulsar's X-ray flux is entirely powered by disk 
accretion onto the neutron star) is ruled out, as it greatly 
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Figure 1. Infrared images of tlie 4U 0142+61 field. 

a, b, Spitzer/IRAC mid-IR images taken on 17 January 2005 
in tlie 8.0 /xm band (a; 75 min exposure) and tlie 4.5 pm 
band (b; 77 min exposure). An 0.15 arcsec astrometric so- 
lution was derived using 14 unblended field stars from the 
2MASS IR point source catalogue. A candidate pulsar coun- 
terpart (indicated by a 4 arcsec radius circle) is detected in 
botli IRAC images at a position of R.A. = 0l''46™22.40^ 
Decl. = -|-61°45'02. 9" (equinox J2000.0), with a la confidence 
radius of 0.7 arcsec. This is consistent with the known pul- 
sar position,^'* lying within 0.3 arcsec. Based on the source 
number density in the 4.5 /im image, the chance coinci- 
dence probability is only 0.2 percent. The source fluxes were 
51.9 ± 5.2 pjy (8.0 /xm) and 36.3 ± 3.6 pJy (4.5 /um), where 
1 ^Jy= 10"^® erg cm"^ s"^ Hz"^. c, A deep Ks-hand (2.2 pm) 
near-lR comparison image of the central part of the IRAC field 
taken on 1 November 2001 from the Keck I telescope in Hawaii. 
(The corresponding subfield is indicated by a box in a and b.) 
Besides the pulsar itself (Ks « 20), there are no other objects 
detected within 4.0 arcsec of the source position down to a 
magnitude limit of Ks ~ 22. The scale and orientation of b, c 
are the same as shown in a. 
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Figure 2. Infrared colour-colour diagram 

for 4U 0142-1-61 (diamond) and 33 field stars (small 
crosses) from the 2MASS point source catalogue. We 

converted the 8 /im fluxes to magnitudes (mg) using a zero 
point of 63.1 Jy, and used published values of the near-IR mag- 
nitudes for the pulsar. ^^ {Ks = 19.9 is an average value for the 
pulsar's variable J^s-band flux.^^) The field stars have colours 
consistent with those of reddened {Av ~ 3.5) main-sequence 
stars (solid line) or red giants (dashed lines). An example 
Av = 3.5 reddening vector is also plotted. The object at the 
position of 4U 0142-1-61 has highly unusual colours (indepen- 
dent of reddening) and a small chance coincidence probability, 
and we conclude that it is the pulsar counterpart. 

overpredicts the optical flux.^^ Instead, we consider a 
passive disk illuminated by the X-ray pulsar. At these 
low temperatures, the disk's continuum emission would 
be entirely from dust. Indeed, we note that our observed 
IR spectrum is quite similar to those of dusty protostellar 
disks, ^ in which the disk is passively illuminated by the 
finite-radius central source. In our case, the neutron star 
is effectively an X-ray point source, but assuming there is 
enough gas to provide pressure support for the dust, we 
would expect any surrounding disk to be slightly flared 
(with corresponding disk thickness h ex r®' '' at disk radius 
r; ref. 18) and thus subject to X-ray irradiation. 

As additional support for an X-ray-heated disk 
model for the IR flux, we note that four of the eight 
known AXPs have known near-IR counterparts (includ- 
ing 4U 0142-1-61), and all of these have similar IR/X-ray 
flux ratios of order 10"'* (ref. 19). In at least one source, 
there is even a detected correlation between the near-IR 
and X-ray fluxes.^" This may indicate that the IR emis- 
sion from all AXPs arises in an X-ray-heated debris disk. 

We fit the power-law-subtracted IR spectrum from 
4U 0142-1-61 with an X-ray-heated disk model^* with 5 
parameters: source distance d, disk inclination angle i, 
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Figure 3. Optical/infrared spectral energy distribu- 
tion of 4U 0142+61. The vertical axes are both scaled by 
frequency v. The left axis shows the //-scaled flux per unit fre- 
quency, vFu\ the right axis shows the //-scaled luminosity per 
unit frequency vLv for an assumed distance of d=3.9 kpc. The 
Ks point is an average value for the pulsar's variable emission 
in that band.^^ The triangles indicate the observed optical/IR 
flux while the squares indicate the dereddened flux assuming 
Av = 3.5. The dereddened optical VRI flux is well flt by a 
Fi, (X p^'^ power-law component (short-dashed line) up to the 
B-band (0.4 /im) break, and is relatively constant in contrast 
to the near-IR variability (0.5 magnitudes). This power-law is 
presumably magnetospheric in origin, an assertion supported 
by the optical pulsations that cannot be due to X-ray repro- 
cessing.^^ The power- law-subtracted IR spectrum is poorly 
flt by a simple blackbody model (long-dashed curve), likely 
excluding the presence of large planets or a very-low-mass 
stellar companion. The best-flt passive debris disk model, 
including irradiation by the central X-ray source, is plotted 
as the dot-dashed curve. The sum (solid curve) of the disk 
and power-law components flts the optical/IR data (below the 
break) well. An accretion disk (including irradiation by accre- 
tion-powered X-rays from the pulsar) can be ruled out because 
it signiflcantly overpredicts the optical flux.^'* 

inner radius rjn, outer radius Tout, and X-ray albedo ?7d. 
This model predicts that the disk's effective temperature 
at radius r is 



T{r) ~5,030K (1 - ryd)^/^ 



3.9 kpc j 



r 
i?0 



-3/7 



(1) 

where an unabsorbed 0.5-10 keV X-ray flux of 4.8 x 10 ^" 
erg cm^^ s^^ from the pulsar (ref. 12) is assumed to il- 
luminate the disk. The total disk flux emitted at a given 
frequency is found by integrating the thermal emission 
over the entire projected (cos«) disk surface between r^n 
and Tout- For our flts, we flxed d = 3.9 kpc and cosi = 0.5 



and assumed r/d > 0.9, as favoured empirically by obser- 
vations of X-ray-heated neutron star disks in low-mass X- 
ray binaries. ^^ We obtained a good fit with Tin = 2.9 i?©, 
rout = 9.7 Rq, and 7]^ = 0.97. 

The inner disk temperature is well constrained to be 
T(rin) 2± 1200 K by the shape of the near-IR spectrum. 
This temperature is comparable to the sublimation tem- 
perature of dust (1000-2000 K, depending upon grain 
composition), suggesting that Tin may be set by X-ray de- 
struction of dust. Disk gas could persist inside rjn without 
any continuum emission at these temperatures, cooling in- 
stead by line emission. Alternatively, rin might also have 
been set by interactions with the pulsar magnetosphere at 
some time in the past, which is consistent with the fact 
that the pulsar's light cylinder radius (rLc = O.GIRq) lies 
well inside rjn . The flt value of rjn is sensitive to our choice 
of rjd] for rjd — 0.5, the disk is further from the pulsar with 
Tin — 19 i?0, although the flt quality is poorer. 

By contrast, the uncertainty in Tout is dominated by 
our ignorance of the spectrum at wavelengths longer than 
8.0 fim, since the location of rout is inferred from the 
frequency below which the spectrum turns over to the 
Rayleigh- Jeans {F^ ex ly"^) regime. (Here Fi, is the flux 
per unit frequency at frequency v.) As this is poorly con- 
strained by our data, our fltted value for rout should be 
regarded as a lower limit. Observations in the far-IR or 
millimetre (MM) bands will be required to determine rout 
precisely. 

Another beneflt of far-IR/MM observations is that 
dust emission in these bands is typically optically thin, 
allowing for a direct measure of the disk mass.^ In the 
absence of far-IR data, we can at least assume that the 
flux at 1 mm does not exceed our observed 8.0-/im flux as 
the spectrum appears to be falling at wavelengths longer 
than 4.5 /im. This sets an upper limit on the total disk 
mass (both dust and gas) of^ 



Afd < 3 X 10"^ Me 
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50/xJy 
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3.9 kpc 
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0.01 cm"^ g 



(2) 



where Mq is the solar mass, i^MM and kmm are the flux 
and opacity at 1 mm, and we have assumed a dust-to- 
gas ratio of 1% in the disk. (Note that our assumptions 
for Kmm and the dust-to-gas ratio do not account for 
the high metallicity and unusual composition possible for 
supernova-processed material.) In protostellar disks, i^MM 
is generally around two orders of magnitude smaller than 
the mid-IR flux,^ and this is consistent with an extrap- 
olation of our best-fit disk model for 4U 0142-1-61. This 
suggests a disk mass of order Md ~ 10~^ M©, compara- 
ble to the total mass of the Earth-mass planets detected 
around the old millisecond radio pulsar PSR B1257+12 
(refs. 22, 23). 

In order to consider the lifetime of the disk Td, we first 
need an estimate of the disk's mass loss rate. For an upper 
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limit, wc can interpret the pulsar's spin period derivative 
(P = 2 X 10^^^) in terms of magnetic 'propeller' torques^^ 
to argue that (— Md) ;$ 10^^^ Mq yr^^. (This presumes 
that the disk gas penetrates inside rin and into the pul- 
sar light cylinder to interact with the magnetosphere at 
the corotation radius Tco — O.Oli?©.) The resulting disk 
lifetime t^ = Md/|Md| /^ 10^ yr is much longer than the 
pulsar's spin-down age Tp = P/2P « 10'"' yr, consistent 
with a supernova fallback origin. This also suggests that 
the disk did not start out with considerably more mass 
than it currently contains. 

Our detection of a fallback disk is, to our knowledge, 
the first direct evidence for supernova fallback in any con- 
text, and it thus bolsters the notion of fallback in a vari- 
ety of other circumstances. For example, direct fallback 
(without a disk) in the seconds to hours following core col- 
lapse has been proposed as a mechanism for black hole for- 
mation following a supernova explosion,^ and may explain 
the absence of a detectable pulsar in the nearby supernova 
SN 1987A. However, if a fallback disk forms around a new 
neutron star, there are several mechanisms that can limit 
the disk lifetime, ^^'^'' including accretion, magnetic pro- 
peller expulsion, and ablation by a pulsar wind, with the 
details depending upon the spin rate and magnetic mo- 
ment of the neutron star. In addition, giant X-ray flares 
such as those seen in some magnetars^^ may also severely 
limit the lifetime of fallback disks. Our discovery implies 
that it is possible for a disk to survive as long as 10'"' yr, 
although 4U 0142+61 may be an exceptional case. 

If fallback disks are common around all young neu- 
tron stars, then pulsar spindown may be caused by disk- 
magnetosphere interactions in addition to magnetic dipole 
torques, rendering traditional estimates of pulsar ages and 
magnetic moments unrealiable.^'^ On the other hand, if 
fallback disks are peculiar to AXPs, this may point to 
a distinct formation mechanism for magnetars. It may 
also explain the narrow clustering of slow spin periods 
for AXPs,^°'^^ despite the need for magnetar activity to 
account for their X-ray flux. However, it is difficult to 
discriminate among these possibilities since there are few 
constraints on the presence of disks around young pul- 
sars. Previous searches, motivated by the discovery of a 
pulsar planetary system, ^^ concentrated on middle-aged 
and old radio pulsars with no success. ^^ Our detection 
in 4U 0142+61 is probably due to a combination of the 
X-ray heating of the disk and the relative youth of the 
pulsar (compared to Td). We note that energetic parti- 
cles in a radio pulsar wind may be as effective as X-rays 
at heating a disk, so that both X-ray and radio pulsars 
are promising targets, and additional debris disks may 
yet be detected in a suitably selected sample. However, 
the strong non-thermal emission from a young pulsar or 
supernova remnant may mask a disk signature in some 
cases. 

By analogy with protoplanetary dust disks around 
ordinary young stars,® the presence of a debris disk 



around a neutron star naturally raises the possibility 
of planet formation. A debris disk is the presumed 
origin^®'^"^ of the system of Earth-mass planets already 
known to exist around the old (^ 10^ yr) millisecond pul- 
sar PSR B1257+12, which was detected through pulsar 
timing residuals. ^^'^^ However, millisecond pulsars are 
generally believed to have been spun up through sustained 
disk accretion, and the formation and survival of a de- 
bris disk and planetary system in this context are dif- 
ficult to explain.'^'' Our discovery does not directly ad- 
dress this system, but it does demonstrate that a debris 
disk around a neutron star is indeed possible. Despite 
this, it is not clear whether the radiation-rich environ- 
ment around 4U 0142+61 is suitable for the formation 
and survival of planets.^'' On the other hand, since planet 
formation probably requires ^ 10® yr after the disk's for- 
mation,'^'' magnetar activity may well have subsided by 
that time.^ In any case, searches for planets around young 
pulsars are likely to be difficult. Pulsar timing searches 
for planets will have poor sensitivity in AXPs owing to 
their slow spin periods, while timing noise in the faster- 
spinning young radio pulsars will limit sensitivity to plan- 
ets in these systems. 
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